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The recent discrepancy between proton charge radius measurements extracted from electron- 
proton versus muon-proton systems is suggestive of a new force that differentiates between lepton 
species. We identify a class of models with gauged right-handed muon number, which contains new 
vector and scalar force carriers at the ~100 MeV scale or lighter, that is consistent with observations. 
Such forces would lead to an enhancement by several orders-of-magnitude of the parity-violating 
asymmetries in the scattering of low-energy muons on nuclei. The relatively large size of such 
asymmetries, O(10 -4 ), opens up the possibility for new tests of parity violation in neutral currents 
with existing low-energy muon beams. 



Introduction. There has been much interest as of late 
in the possibility of new gauge forces existing in the MeV- 
GeV scale, stimulated in part by the prospect of a light 
mediator between dark matter and the standard model 
(SM) (see, e.g., [l|). While many models of this type 
can be explored, a great deal of attention has been given 
to a new /7(1_) gauge boson V kinetically mixed with 
hypercharge [2] . At low energies V appears as a massive 
copy of the ordinary photon, 
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where k is the mixing angle parameter. The conservation 
of the electromagnetic current and the absence of any 
intrinsic parity/flavor/CP violation in the interaction of 
V with the SM fermions can hide this force from very 
powerful symmetry tests. The model ([1]), while perhaps 
the simplest, is not the unique possibility for new gauge 
interactions below the weak scale [3] . 

While the astroparticle physics incentives are rather 
speculative, an additional motivation for a new light 
gauge boson is provided by terrestrial experiments. 
Among several discrepant low-energy measurements, the 
recent determination of the proton charge radius using 
muonic hydrogen [i| and the long-standing measurement 
of the anomalous magnetic moment of the muon [H[ may 
be manifestations of a new sub-GeV scale force carrier 
that couples preferentially to muons. In this letter we 
argue that if such discrepancies are caused by a new 
muon-specific gauge force, one should expect parity non- 
conservation (PNC) in the scattering of muons on nuclei 
far above the SM level. We point out the feasibility of 
a dedicated search for the PNC asymmetry, enhanced to 
O(10 -4 ) level, with existing low-energy muon beams. 

With our present understanding of the strong inter- 
actions, the charge radius of the proton r p cannot be 
computed from first principles but instead must be ex- 
tracted from experiment. The comparison of r p values 
obtained using different experimental methods provides 
a consistency check of QED theory and constrains a va- 
riety of new physics scenarios. Currently, there are three 



competitive ways of determining r p : 1) high-precision 
measurements of the atomic levels in hydrogen and deu- 
terium, 2) direct electron-proton scattering experiments, 
and 3) the measurement of the Lamb shift in muonic hy- 
drogen. The most precise determinations currently read 



r pA = 0.8768(69) fm 
r p , 2 = 0.879(8) fm 
r p , 3 = 0.84184(67) fm 



atomic H, D, (2) 
e — p scattering, (3) 
muonic H. (4) 



The r p values obtained from e — p systems are consistent 
with each other and significantly differ from the r p value 
extracted from the muonic hydrogen Lamb shift, 



r P ,i ^ r P)2 > r Pj3 , 
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The difference between rj>i @ and r Pi 3 [1] is 5a while the 
difference between r p ^ [7| and r Pi 3 is 4.6cr (for an up-to- 
date theoretical analysis see Ref. [§[). Part of this dis- 
crepancy may be related to the model dependence of the 
proton form factor used in various extractions of r p 0], 
and it is conceivable that further scrutiny of SM predic- 
tions can close this gap [loj]. At the moment, however, 
this discrepancy stands and has stimulated investigations 
of new interactions that could potentially be responsible 
for the difference [IlT[l3j ] . The difficulties associated with 
such an enterprise stem from the fact that the difference 
([5|) requires the strength of the new interactions to be 
on the order of O(10 4 Gf), which is impossible to attain 
without new light states below 1 GeV. 

It is easy to see that the kinetically mixed vector ([T]) 
cannot explain the observed pattern. In the presence of 
V-exchange, the inferred r p would actually depend on 
the effective momentum transfer |g| involved [14J. For 
r p,i ( r p,3)> this corresponds to the inverse Bohr radius 
~ am e (am^), while for r p ^ the momentum transfer is 
much larger. The effect of the extra attraction generated 
by V will be interpreted as the largest negative correction 
to the charge radius for the experiment that involves the 
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smallest \q\. Therefore, a kinetically mixed vector pre- 
dicts r Pj i < r p .3 < r p .2, which is not consistent with the 
observed pattern ([5]) . One can easily show that the inclu- 
sion of several kinetically mixed vectors does not change 
this pattern. Another logical possibility is a repulsive 
Yukawa force between protons and muons/electrons, e.g., 
as may occur if there is a new force with gauged baryon 
number and kinetic mixing with photons. However, in 
this case the natural pattern will be r p ^ < r p ^ < r p> i, 
which again disagrees with ([5]). 



In Refs. [12J, 1 1 31 ) a purely phenomenological approach 
to explain ([5]) was taken, in which dimension 6 operators 
(p,j a fi) (p'jaP) or (fill) (pp) are mediated by the exchange 
of a new light vector or scalar particle. Scalar mediators 
of this type are reminiscent of a very light Higgs boson 
and will face stringent constraints from rare meson decays 
and neutron-nucleus scattering 15(. Vector mediators 



are more promising, but in order to be integrated with 
the rest of the SM, the following conditions must be met: 

i. The interactions must be formulated in terms of SM 
fermion representations, 

ii. No new interactions stronger than Gf can exist be- 
tween neutrinos and nucleons or electrons, 

Hi. No new electrically charged elementary particles with 
masses below 100 GeV can exist, 

iv. The model must have the possibility of a UV comple- 
tion at or above the weak scale, 

v. The model must be consistent with a variety of tests 
from QED and particle physics in the MeV energy range. 

The second condition comes from the wealth of data on 
neutrino scattering in the E ~ 10 MeV energy range and 
neutrino oscillations and is emphasized here because it 
serves as a powerful model discriminator. Indeed, the 
interaction of a new particle V with the lepton vector 
current may be viewed as a subset of the interaction with 
left- and right-handed SM fermion currents, 

V a lj a l C V a (cxLrf a L + C2Br( a R), c\ ^ -c 2 . (6) 

The left-handed fermion doublet L includes a neutrino 
field, so the requirement ii. is equivalent to c± = 0. This 
forces V to couple to the pure right-handed fermion cur- 
rent. The absence of large neutral right-handed currents 
for electrons follows from PNC tests in the electron sec- 
tor, and we therefore conclude that the most promising 
coupling of a vector particle that can explain ([5]) is to the 
right-handed muon. 

Models with gauged \ir. We now focus on the class 
of models based on a new U(1)r gauge symmetry with 
quantized fiji number. The Lagrangian is 
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Here V is the U(1)r gauge boson, (f> is a new Higgs field, 
neutral under the SM gauge group and charged under 
U(1)r, that condenses (0) = vr/^/2, D = d + ignQsY + 



ieQEAiA, and k is the mixing angle parameter. The mass 
term for the muon is necessarily a higher-dimensional 
operator involving both </> and the SM Higgs field Hg M 
generated at a high scale A, 
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with vh/(V2A) entering as an effective SM-like Yukawa 
coupling for the muon. As we shall see below, the range 
for vr suggested by the charge radius phenomenology is 
fully consistent with A being at the weak scale. There- 
fore, we are not concerned with building an explicit model 
that provides a UV completion to C m . The physical ex- 
citation of is a new muon-specific Higgs scalar S in the 
mass range ms < vr. 

The model ([7]) suffers from gauge anomalies involving 
the photon and V. It is possible to restore gauge invari- 
ance by introducing dynamical scalar 'gauge' degrees of 
freedom. The price for maintaining gauge invariance is 
that the theory becomes nonrcnormalizable, with a UV 
cutoff Ajjv above which calculability of the theory is lost. 
The estimate for Kjjv may be obtained, e.g., from the ra- 
diative three loop vector self energy diagram [l6| : 
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We observe that vectors in the range my ~ 10 — 100 
MeV with couplings qr ~ 0.01 — 0.1 are consistent with 
a UV cutoff well above the TeV scale. There are of course 
examples of perturbative cancellations of the anomalies, 
such as quantized [1r + sr — cr. Such a scenario faces 
severe constraints from quark flavor physics, cc resonance 
decays, and parity-violating tests involving nucleons and 
appears to be thoroughly excluded. Therefore, we choose 
the model ([7]) as the best candidate to describe new 
muon-specific forces, which is a consistent effective field 
theory valid below Ajjv ©■ 

Phenomenological constraints. From the Lagrangian 
(Jjj we obtain the couplings of the new vector and scalar 
particles to fermions, 
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9A P ' n = 9v = °; 9s = \9R\ m n/ m v, 



where e = \/4ira is the positron charge. With the cou- 
plings (fTU)) . we calculate the corrections to the energy 
levels of the ordinary and muonic hydrogen that will be 
interpreted as corrections to the proton charge radius. 
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where a = (am^mp) (m^ +m p ) is the /iH Bohr radius, 
f(x) = x 4 (l + £)~ 4 , and r\ = Kgn/(2e). The difference 
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must be consistent with the observed 
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Parameter 


romt A rOint ±3 romt \u 


mv 


10 MeV 50 MeV 100 MeV 


m s 


~l r\C\ O A Ti T XT t~\ r\ A A TV T T T Ci A /"\^T 7\ !f T T 

102.84 MeV 90.44 MeV 84.97 MeV 


9R 


0.01 0.05 0.07 


K 


0.0015 0.0075 0.02 


V 


2.5 x 10~ 5 6.2 x 10~ 4 2.3 x 10~ 3 


Vr 


1 GeV 1 GeV 1.4 GeV 



TABLE I: Benchmark points for the model that pass all phe- 
nomenological constraints. 



pattern (|5|) and requires rj to be positive. In the scaling 
regime of amy S> 1 one has 
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2.5 x 10~ 5 
(10 MeV) 2 
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In the same regime, the model predicts that future exper- 
iments with /^He would detect the effective charge radius 
of the helium nucleus shifted down by Ar^ = — 0.06fm 2 . 

Another important constraint comes from the mea- 
surement of g — 2 of the muon, which currently displays a 
~ 2 — 4cr discrepancy with the SM prediction depending 
on the estimate of the hadronic contribution. The in- 
dividual one-loop contributions from V and S are large 



(i.e. enhanced compared to the pure vector case [14 1 



by m^/my) and of opposite sign, so that for the choice 
of parameters (fT2"|) they must cancel. This mutual can- 
cellation must happen at a per-mille level, and should 
be considered as the main phenomenological drawback 
of the model ([7]). Nevertheless there do exist choices of 
parameters that simultaneously account for the charge 
radius data and the muon g — 2 discrepancy. 

Other constraints that must be taken into account are 
the electron g — 2 determination vs. independent mea- 
surements of a 11811 and tests of d — p transitions in 
muonic Si and Mg [l9j , for which no deviations from stan- 
dard QED predictions were found. Table I displays three 
benchmark points for my = 10, 50, 100 MeV for which 
all constraints are satisfied. Vector masses my < 10 MeV 
are excluded by muonic Si,Mg data and tests of a. 

Additional constraints on gauged /j,r theories depend 
on the decay channels of S and V. If no new states 
charged under U(1)r exist below my/2, the gauge boson 
V will decay to e + e~ pairs and thus be subject to tests 
at lepton colliders and fixed target experiments [2(J. In 
particular, a preliminary search for the rare decay mode 
4> — > r/V would disfavor models with k ~ O(10~ 2 ) and 
my above 30 MeV 21] . If new decay channels for V are 
allowed these bounds can be relaxed. Among model in- 
dependent probes, the different couplings of V to muons 
vs. electrons (1101) suggest nonuniversal leptonic branch- 
ings o f J/ ip and other narrow vector resonances. Current 
data 22J is only sensitive to r\ > O(10 -2 ), which does 
not probe the most interesting my <J 100 MeV regime. 
An alternative way to search for ^-exchange is to study 
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FIG. 1: The asymmetry A L r(6>) defined in Eq. for the 
benchmark points labeled A, B, and C in Table |U The solid 
curves are for p — 29 MeV/c and dashed curves for p — 
200 MeV/c. 



the 0(r)) forward-backward asymmetry in e + e~ — > fi + \x~ 
annihilation at medium energy high-luminosity facilities 
with longitudinally polarized beams. Finally, since the 
gauge boson couples preferentially to //r, there is no rea- 
son for charged lepton flavor conservation. One must 
therefore make some assumptions about the underlying 
structure of the Yukawa matrices, implying new flavor- 
related physics in the UV. 

New parity-violating effects. Despite the existence of 
polarized muon sources, no tests of PNC in neutral cur- 
rents involving low-energy muon beams have been per- 
formed. This is because the maximum muon intensity 
corresponds to p — 29 MeV/c, where the parity violating 
asymmetry due to the weak interactions will not exceed 
O(10 -7 ). With the introduction of a new vector force 
coupled to jttfl, the PNC effects are greatly enhanced. 
For the scattering of semi-relativstic muons on a heavy 
nuclear target, the asymmetry is given by 



dcr L - da R 
dai, + dart 
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Q 2 + m 2 v l-f3 2 sin" {9/ 2) 



(13) 

where Q is the momentum transfer of the elastic scatter- 
ing, Q 2 /p 2 = 4sin 2 (6»/2), /3 = \p\/E, and L(R) label the 
incoming muon's helicity. Notice that the same combi- 
nation of couplings 1] governing the correction to r p also 
determines the asymmetry. The PNC asymmetry j4lr(#) 
is presented in Fig. ([1} for the three benchmark points in 
Table U with two reference values of incident muon mo- 
mentum, 29 and 200 MeV. The asymmetry can vary in 
a broad range from 10~ 5 to 10~ 3 , becoming larger for 
Q 2 > m 2 , due to the scaling relation (ji"2"j) . 

We next investigate the feasibility of achieving sta- 
tistical sensitivity to Alr in Eq. (fl3|) in a Rutherford- 
type scattering setup. Since low-energy muons are easily 
stopped, counting rates are maximized with the use of 
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high Z thin foil targets, while the optimal Z should be 
determined from the combined analysis of statistical and 
systematic errors. If, for example, a tungsten (Z = 74) 
foil of d = 0.01 mm thickness is used, the muons will lose 
only ~ 5% of their kinetic energy. Assuming a muon- 
counting detector with full azimuthal coverage and polar 
angle coverage in the range from 60 to 80 degrees where 
the asymmetry is maximized, one obtains the following 
probability for the scattering of a muon at a large angle: 



P = dx N atoms x Volume 1 x CT Rth - 6 x 10" 



(14) 



With this probability, the time required to collect N 
(Alr)~ 2 events is given by 



*|jV~10 s 



N 
PC 



1600 s x 



10 s muons/s 



(15) 



where we have normalized the muon flux to the highest 
modern beam intensities 23[. It is thus apparent that 



the statistical uncertainty will not be a limiting factor in 
detecting parity violating asymmetries of order 10~ 4 . 

Another promising avenue in the search for anomalous 
PNC effects is the study of parity-violating decays of 2s 
states in muonic atoms, in which PNC will manifest in 
the enhanced one-photon rate of 2s decays. To illustrate, 
we assume rj/m v is fixed by the scaling limit (fT2|) and 
compute the 2sx/2-2px/2 mixing in /i 4 He. The results for 
the mixing angle 6, ratio of El to Ml amplitudes for 
the one-photon decay of the 2s!/ 2 state, and rate for the 
one-photon decay are given by 



6~ 3 x 10~ 5 ; 



A{El) 
A{M1) 



60; 



(16) 



r^i, ^ 1.9 x io 3 Hz; r^ ls /r^ ls ~ o.ois. 



The rate of the one-photon decay is only marginally 
smaller than the one-photon quenching rate 24) at a gas 
pressure of 4hPa, at which the /i 4 He Lamb shift exper- 
iment is planned, and the presence of such decays can 
be searched for at different gas pressures (with a modest 
improvement of the 27 background rejection). 

Finally, despite the absence of parity- violating cou- 
plings to the electron at tree level, an e-nucleus parity- 
violating amplitude will still occur at the two loop level. 
Given the accuracy achieved in tests of PNC with elec- 
trons, it is therefore highly desirable to calculate this 
effect, which may lead to an independent constraint on 
gauged /i/j models. 

To conclude, we have argued that the class of models 
with gauged [Ir represents one of the few possibilities in 
which the discrepancy between e— p and fj,— p determina- 
tions of the proton charge radius can be reconciled with 
new physics. Although anomalous, these models consti- 
tute valid effective field theories which can in principle 
be UV-completed at the weak scale. The simultaneous 
explanation of the r p data and muon g — 2 discrepancy 



requires a tight correlation between the scalar and vector 
masses. A striking consequence of this class of models is 
the existence of enhanced PNC effects in the muon sector 
that can be searched for at existing muon beam facilities. 
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